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PETER C. GIBSON 

Abstract. We analyze the correspondence between finite sequences of finitely sup- 
ported probability distributions and finite-dimensional, real, symmetric, tridiagonal 
matrices. In particular, wc give an intrinsic description of the topology induced on 
sequences of distributions by the usual Euclidean structure on matrices. Our results 
provide an analytical tool with which to study ensembles of tridiagonal matrices, 
important in certain inverse problems and integrable systems. As an application, 
we prove that the Euler characteristic of any generic isospectral set of symmetric, 
tridiagonal matrices is a tangent number. 
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1. Introduction 

In this paper we analyze the correspondence between finite sequences of finitely 
supported probability distributions and finite-dimensional, real, symmetric, tridiag- 
onal matrices. In particular, we give an intrinsic description of the topology induced 

Research at MSRI was supported in part by NSF grant DMS-9810361. The author's research is 
supported also by NSERC Postdoctoral Fellowship 231108-2000. 
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on sequences of distributions by the usual Euclidean structure on matrices. Our 
results provide an analytical tool with which to study ensembles of tridiagonal ma- 
trices, which are important in a variety of contexts, including inverse problems 
and integrable systems [|T|. 

The paper is divided into two parts. In the first part. Section]^, we study the map / 
that takes the spectral distribution of a finite Jacobi matrix to that matrix. The open 
d-simplex 5" naturally represents probability distributions whose support is a fixed 
{d + l)-element set A of reals, and it is well-known that / maps S homeomorphically 
onto the manifold M of Jacobi matrices having spectrum A. However, viewed as a 
map on the closed simplex S, f has points of indeterminacy on the boundary dS. 
Our main result is the explicit construction of a blow-up {B,7r) of S, centred at dS, 
that eliminates the points of indeterminacy of /, establishing a homeomorphism h 
from B onto the closure of M. 

The original motivation was to better understand geometrically solutions to certain 
ill-posed inverse problems concerning Jacobi matrices, specifically p|. Theorem 6.4] 
and [p. Theorem 2.3]. The essential question. Problem 1 on p.|^, is to describe intrin- 
sically the topology induced on sequences of distributions by their correspondence 
with direct sums of Jacobi matrices. A complete answer to the question is implicit in 
the blow-up {B, vr) itself. As far as we know, these results are new; their application 
to inverse problems will appear in a future article. 

The second part of the paper. Section]^, consists of an application of the first part to 
isospectral sets of real, symmetric, tridiagonal matrices. It turns out that B, although 
not convex, has a facial structure like that of a convex polytope. Based on this, we 
use the homeomorphism h : B ^ M to construct a family of polyhedral complexes Pd 
that encode the topology of the isospectral sets in question. This enables us to prove 
that the Euler characteristic of a generic isospectral set of symmetric, tridiagonal 
matrices is a tangent number. From the point of view of discrete geometry, the 
complexes Pa are remarkable in their own right. They have a simple description, yet 
exhibit exotic combinatorial properties. 

The topology of isospectral sets of tridiagonal matrices has been considered inde- 



pendently, and much earlier, by Tomei, Nanda, Deift, Flaschka and others [|T0[, 0, 
. But whereas these authors make essential use of the Toda flow, our methods are 
more purely combinatorial. In place of a dynamical system, we exploit the above- 
mentioned blow-up construction, which is completely explicit. Tomei also computed 
Euler characteristics, but does not seem to have noticed that hyperbolic tangent is a 
generating function for them. In any case, it is of interest to try to understand what 
connections exist between the two apparently different approaches. 

Acknowledgement. Thanks both to Richard Cushman and Jiirgen Bokowski for 
many helpful conversations at various points in the preparation of this article. Thanks 
also to Gunther Uhlmann and MSRI for the stimulating research environment of the 
MSRI semester in inverse problems. Fall 2001. 
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2. Sequences of distributions, direct sums of Jacobi matrices 



Probability distributions proved useful in for parametrizing solution sets of 
certain classes of inverse problems concerning finite Jacobi matrices, based on the 
well-known correspondence between distributions, orthogonal polynomials and Ja- 
cobi matrices — see . In the present section, we study an extension of the standard 
correspondence, considering sequences of distributions and direct sums of Jacobi ma- 
trices in place of, respectively, distributions and Jacobi matrices. As will be explained 
in a future article, this extension gives a clearer geometric interpretation of the above- 
mentioned solution sets. 



In Sections and p.2| we assemble the standard facts, without proof, and establish 
the basic notation to be used throughout Section A precise formulation of the main 
problem is given at the end of Section With the exception of the formula cited 
in Proposition |l^ of Section |2.6| , the material in Sections |2.3| - |2.8| is new and so we 
give detailed proofs of all results. 



2.1. Notation, basic facts, and definitions; Problem 1. We are interested in the 
correspondence between probability distributions on the real line and Jacobi matrices, 
restricted to the finite-dimensional setting. With respect to this correspondence, 
sequences of orthogonal polynomials play a useful auxiliary role. 

Let JT" denote the set of finite-dimensional Jacobi matrices, that is, matrices of the 
form 
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for some d > 1, where each a„ G R, and each a„ that has even index n is strictly 
positive. We label the entries of J as in (|I|) in order to suggest the map 

(2) J (ai,a2, . . . ,02d-i), 

whereby d x d Jacobi matrices correspond to points in R^'^"^. 

Let V denote the set of finitely supported probability distributions on the real line. 
Elements of V are most simply represented by combinations of translates of the Dirac 
delta function, of the form 



(3) 



da{\) = ^m„(5(A - A„), 



n=l 



where each A„ G R, each w„ > 0, and Yln= 



= 1. For technical reasons we shall 
adopt a slightly more cumbersome, equivalent representation, whereby the weights 
Wn in (PD function as homogeneous, rather than absolute, coordinates. More precisely. 



let V denote the set of distributions of the form 



where each \n,uJn G R and 
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each WjnWn > 0. For da, d(3 G V , write 

da ~ df3 

if da = rdjS for some non-zero r G R. Then an element of V is represented by the 
equivalence class modulo ~ of a distribution in V, consisting of its scalar multiples. In 
symbols, V = V / ~. Having established the precise meaning we will be less formal in 
practise; given da G V, we will simply write da G V, letting the representative stand 
in for its class. The point of this notation is that we are free to rescale probability 
distributions without changing them. 

Let MOP denote the set of finite sequences of monic, orthogonal polynomials. 
These are sequences of the form 

{Po,Pu ■ ■ ■ ,Pd), 

where d > 1 and each Pn{^) G R[A] is monic of degree n, and for which there exists 
a probability distribution da on R such that 

/oo 
plda > for each n, 
-oo 

/oo 
PmPn da = if m 7^ n. 
-oo 

Each J E J corresponds to a unique da G P, its spectral distribution function, as 
follows. Orthogonally diagonalize J as 

J = OAO^, 

where 00^ = I and A = diag(Ai, . . . , \d), and set 

d 

(4) daj{X) = Y,{0{l,n)f6{X-Xr.), 

n=l 

where 0{i,j) denotes the (z, j)-entry of the matrix O. This is the (normalized rep- 
resentor of the) spectral distribution of J in the sense that for any analytic function 

/, 

/oo 
fdaj. 
-oo 

An equivalent characterization is the formula 

(5) r^M^. 

J-oo Z — A 

The map J ^ daj is a bijection from J' onto V. Our main focus in the present 
paper is its inverse 

Here we define / explicitly, along with an auxiliary map 

* : P ^ MOP. 
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Definition 1. Let da be given, and write d = \suppda\. Define 



(6) 



f{da) 
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■ ■ ■ ,pd), 

and {po, 



(7) ^{da) = {po,Pi, 

where the sequences (ai,... ,a2d-i) 
follows. 

Initial step For the purpose of notation set ao = 0. Set 

(1, ^)da 



. ,pd) are constructed recursively as 



Po(A) 



ai 



Pi{X) = X-ai. 



{lA)d^' 

Continuing step Given a2n-2, CL2n-i ond pn-i,Pn, for n in the range 1 < n < d — 1, 



set 



(9) 
(10) 



02n 

02n+l 
Pn+l(A) 



{Pn,Pn) da 

da 

{Pn, Xpn)da 
da 

(A - 02n+l)Pn{X) - oj^Pn^liX). 



The sequence of polynomials "^{da) has a direct relationship to the matrix J = 
f{da), namely, p„ = \E'„((ia) is the characteristic polynomial of the leading n x 
n submatrix of J; in particular, pd is the characteristic polynomial of J. (This 
fact will be needed in technical arguments later on.) Also, the truncated sequence 
[po, . . . ,Pd-i) is orthogonal with respect to the starting distribution da, but 

/■oo 

p'^da = 

is degenerate. Indeed supp da is the set of roots of pa- If one extends J to a (rf + 
1) X {d+ 1) matrix J' by choosing arbitrary entries a2d, 02d+i, then the full sequence 
{po, . . . ,pd) is orthogonal with respect to the spectral distribution dajr of J'. We 
index the components /„ of / according to the sequence (ai, . . . , a2d-i), so that 

On = fn{da). 

Let us consider a fixed d = | supp da |, and write 

d 

da{X) = ^m„(5(A - A„). 

n=l 

Observe that, according to Definition |l|, for 1 < n < d — 1, 

• and /2n+i are rational functions in the variables wi, . . . , Wd, Ai, . . . , A^; 

• /!„ and /2n+i are homogeneous functions of degree in the variables wi, . . . , Wd- 
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Now, / itself is not rational, because of the square root (^). We have chosen the 
positive square root, and in some sense this choice is artificial. It is more natural 
to allow both square roots, and thus to make / a multi-valued function (precisely, a 
2''~^-valued function for fixed d = |supp(ia|). This point of view is developed fully 
in Section |; for now we leave / as single- valued. 

Having determined the nature of the functions /„ in terms of the weights w„, let us 
briefly consider the reverse question, restricted to to matrices J having a prescribed 
set {Ai, . . . , Xd} of eigenvalues. To that end, let 

d 

(11) y^^WnSjX - An) = da J 

n=l 

denote the normalized spectral distribution of J, as defined in (|). Then the nor- 
malized weights Wn are polynomials in the entries ai, . . . ,a2d-i of J, a fact easily 
inferred from (^, as follows. Let J^^ denote the submatrix of J obtained by deleting 
the first row and column, and let 

d 

PdW = l[^\-Xn) 

n=l 

be the characteristic polynomial of J. Then the left-hand side of (^) has the form 

det{zl-j'')/pd{z), 

while the right-hand side of (|^), expressed in terms of (|TT]), is the partial fraction 
decomposition 

d 



Comparing these two expressions yields 

Wn = det(A„/ - j")/|)'rf(A„), 

which, for fixed Ai, . . . , A^, is a polynomial in the entries of J. 

The set of d x d symmetric, tridiagonal matrices is a subspace of R*^^, isomorphic 
(indeed isometric) to R^*^"^ by the map (^. With this isometry in mind, we will 
usually write J G R^"^"^ when J is a. d x d symmetric, tridiagonal matrix. Now, 
because J'nR^'^"^ is not closed, it is useful to consider the extended class of matrices 

j= |Jci(jrnR2'^-^). 

JT" consists of direct sums 

J = Jl © ■ ■ ■ © Jr 

of Jacobi matrices Ji, . . . , Jr & J^, where r > 1. By analogy we extend V to the class 
V consisting of all sequences 

dai © ■ ■ ■ © dar 
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of distributions dai, . . . ,dar G V, where r > 1. (We could alternatively write 
{dai, . . . , dttj.) for a sequence of distributions; the symbol © just serves to emphasize 
the connection with matrices.) The map 

f:V^J 

extends to a map 

in the obvious way, namely, by defining 

J{dai © ■ ■ ■ © dar) = f{dai) © ■ ■ ■ © f{dar). 

Note that, by its very definition, / is a bijection. 

For each d > 1, J' Ci R^*^"^ inherits the Euclidean structure of R^"^"^. But it is 
unclear a priori what it should mean for two points, 

rftti © ■ ■ ■ © rfa^, d(3i ® ■ ■ ■ ® d(3s, 

in V to be "close" . Thus, roughly speaking, a basic problem is to describe intrinsically 
the topology which D acquires from J' via /. We make this more specific as follows. 

Fix a (i-element set A C R. Denote by J' the set of d x d matrices J & J having 
spectrum A, 

= [J \ J ^ Rad-i and spec J = A } . 

Write to denote the corresponding set of sequences of distributions, that is, se- 
quences 

dai® ■ ■ ■ ® dar eV 

such that 1 < r < (i, 

r r 

(12) |supp dan\ = d and \^ supp dan = A. 

n=l n=l 

With this notation, the restriction of /, 
is a bijection. 

Problem 1 Describe intrinsically the topology which J'^ induces on via /. That 
is, characterize convergent sequences of points in V , and show how to compute limits 
of such sequences without evaluating /. 

We illustrate Problem 1 with a simple example. Let dao, dai G V denote the 
distributions 



dao^X) = 5(A), 

dai{X) = (5(A-l) + (5(A-2), 
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and consider the hmit 



L = f 




hm f{dao + t dai) 



Given that the hmit exists, what is LI A naive guess might be L = da^ © dai, but 
this is wrong, as can be verified numericahy by coding the recursive definition of / 
(in, say, Matlab) and then computing f{daQ + tdat) for decreasing values of t. In 
fact, L = dao © da'i, where 



In other words, for A = {0, 1, 2}, the curve dPt = dao+tdai (or equivalently, sequence, 
if we set t = l/ra, n G Z+) converges in V , with respect to the topology induced by /, 
to dao ffi dai. And although dPt £ 1^ for every t > 0, its limiting value as t — >• 0"*" 
is not in V. Rather, df3t "splits" in the limit into a pair of smaller distributions. 
Problem 1 calls, in part, for a full description of this phenomenon. 

2.2. Preliminary observations, further notation. Let A C R be a d-element 
set, for some d > 1, to be held fixed throughout the present section. As a first 
step toward solving Problem 1, we note that V is a disjoint union of open discs of 
various dimensions. The essence of Problem 1 then consists in describing how these 
component discs fit together. 

The notion of an ordered partition is useful for present considerations. Given a set 
S, a sequence 



of sets is an ordered partition of S if: each S'„ is non-empty; Sm fl S'„ = whenever 
m ^ n; and IJn=i = S. Denote the set of ordered partitions of 5* by Vs- The 
condition (|12D for dai © ■ ■ ■ © dar to belong to V is precisely that 



be an ordered partition of A. And for each ordered partition (5*1, . . . , Sr) G V\, there 
corresponds a piece of V , 



da[{\) 



5(A- 1) +45(A-2) 
X'^dai{X). 




(supp dai, . . . , supp dar) 



(13) {dai (B ■ ■ ■ (B dar ^ T) | (supp (iai, supp (ia^) = (5*1, . . 
which is isomorphic to an open disc of dimension 



,Sr)}, 



r 



n 



1 



d — r. 



n=l 




A„) © ■ ■ ■ © WnSiX 
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the weights m„ of which serve as coordinates on (p!3|) . Indeed the notation is simplest 
if we fix an indexing of A, 

A = {Ai,... ,A4, 

and then work with partitions of 

[rf] = {!,... ,4 

rather than of A. We set up notation for our coordinate space as follows. 
Given a non-empty set 5 C [d], we write for the projective space 

= {w:S supp w ^ 0}/ ~, 

where w ^ w' if w = rw' for some non-zero r G R. Thus consists of sequences 
indexed by S, distinguished up to proportionality, and is isomorphic to RP'"^'"^. It 
is important to distinguish from whenever S ^ S'. When it is desirable to 
emphasize the domain of w : S* ^ R, we append it as a superscript, writing . As 
usual, we write Wn, or w^, to denote the value of w at n e 5*. We will always use 
homogeneous coordinates on W^, treating w G as a particular map w : S' — > R, 
and letting the equivalence relation ~ be understood. We write Wf for the open 
subset of defined by 

= {w: S I WmWn > Vm, n e S} . 

Referring to ([T3|), let {S[, . . . , S'^.) G V[d] be the partition corresponding to (5*1, . . . , Sr) G 
Va, so that for each n, 

S'^ = {te[d] I A, g5„}. 
Then the weights Wn in (0) belong to the coordinate patch 
(15) '^^^ = 1^+^ © ■ ■ ■ © W+^, 

and we may identify (|13|) with W_^_ ^' "' . We denote an element 

using the abbreviated notation 

^{Sl,...,Sr) ^ ^ . . . ^ ^S,_ 

(S' S') 

Now, WP'-' ^' is a direct product of open subsets of finite-dimensional projective 
space, and carries its own topology. By the natural topology of (|T3|), we mean the 

topology of this direct product. We identify V with the union of coordinate patches 
([TSl) , indexed by partitions of [rf], 

(Si,...,5.)GP[dj 
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and use the obvious notation for the function / : — > j'^ . Namely, for G W^, 
f{w^) denotes 



\n&S J 



and, given a point w^^^''"'^''^ G pv^j!^^'"' we write 

It is easy to see that the topology on induced by /, restricted to a given piece 

( S S \ ( s s ^ 

agrees with the natural topology. The restriction to W^. '' of each 
component /„ of / is either a rational function without singularities, or the square 
root of such a function, and the components Wn of the normalized representor of 
the of the restriction of / to + ) are polynomials in the matrix entries 

ai, . . . , a2d-i- Therefore the restriction to \Y^^''"'^^^ of / is a homeomorphism. But 
there is no natural topology on V as a whole; we need to analyze / directly to 
determine how its pieces fit together in the induced topology. Note that V has a 
single piece, W^f, of maximum dimension d — 1, corresponding to the trivial ordered 
partition {[d]) G V[d]. We refer to w!l^ as the principal component of . It turns 
out that the induced structure of is determined by sequences, or curves, in the 
principal component which tend to the boundary of in W^'^h In fact it is enough 

just to consider moment curves in w]^ tending to the boundary. Note that w]l^ may 
be thought of as the standard, open, {d — l)-dimensional simplex, as each point in 
has a unique representative w satisfying X]n=i ""^n = 1, with each w„ > 0. 

2.3. A blow-up of the standard simplex. Given a (i-element set 

A = {Ai,... ,Ad}cR, 
we construct a pair {B^, vr), consisting of a semi-algebraic set and a bijective map 

The set B^ is defined in a relatively high-dimensional projective space W by quadratic 
equations and inequalities. The W-structure of B^ encodes the topology on V 
induced by / in that the map 

Jon:B^^J^ 

is a homeomorphism. (This latter fact will be proved in Section p.7| .) The pair {B^, n) 
is a blow-up of w^'^'--- '^''^ c. V with centre 

]j{Si,...,Sr) ^ ^j^g boundary of i'-'^-) in yyCSi,... ,5.)^ 

for every ordered partition (5*1, . . . , Sr) G in the sense that there is a coordinate 
projection 7r('5'iv,<S'r) q,^ yy such that: 
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• the restriction of vr to tt coincides with -ji^^^'---^^'-) _ 

In particular, vr) is a blow-up of the standard simplex w!^^'^ C V^. 
We now describe {B^, vr) in detail. is contained in the ambient space 

The dimension of W is easily computed to be {d — 2)2^^"^ + 1. Given w G W, we 
write for its coordinate projection onto the factor , and for the entry of 
corresponding to n & S. Note that, as established earlier, these are homogeneous co- 
ordinates. The equations defining B^ involve the following scalars whose provenance 
we shall explain later, in Section |2]^. For each non-empty S C [d] and each i G [d], 



set 

= l[{>^n-X^y and Cf = l. 

Definition 2. Let B^ denote the semi- algebraic subset o/W consisting of points w 
which satisfy the equations and inequalities 

(16) Ct^wfwf = C^^^wfw^ Vi?C5C[d], i,jeR 

(17) wfw^ > V5C[rf], i,jeS. 



We need some technical preparation in order to define the map tt : B^ V . Here 
it will be useful to consider the collection of subsets of [o?], partially ordered by 
the inclusion relation C; we make use of a correspondence between totally ordered 
sets, or chains, in C[d] having maximum element [rf], and ordered partitions of [d]. 
Let w G S* C [d],n E S, and write P{S, n; w) for the statement 

w^^O and VT C [d],wl ^ ^ |r| < \S\. 

Lemma 3. // P{R, m; w) and P{S, n; w) then either R C S or S C R. 

Proof Set T = R U S and choose i E T such that wf ^ 0. (Such an i exists 
because are homogeneous coordinates for the projective space W^"^.) If z G -R, 
then Cf~^w^wf ^ 0, so by equation (0), C'^^wf'w^ ^ 0, which implies iv^ ^ 0. 
Therefore, since P{R,m]w), \T\ < \R\, forcing T = R, or equivalently, 5" C i?. 
Similarly, if z G S*, then P{S, n; w) and equation (|l^) force i? C S*. I 
It follows from Lemma ^ that for each n G [d] there is a unique S C [d] such that 
P{S, n; w), so that we may define 

'■ [d] 

by setting 

'^w{n) = the unique S C [d] such that P{S,n; w). 

Note that [d] itself is in the range of y^^, since necessarily wf^ ^ 0, for some 
whence P([c?],i;w) automatically holds. Moreover, again by Lemma ^, the range 
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{Ki, ... ,Kr} C of (puj is totally ordered; i.e., for an appropriate relabeling it is 
a chain of the form 

(18) [d\ = Ki D ■ ■ ■ D Kr (H, 

for some r > 1. We are now in a position to define the map vr, but first we state 
another technical fact, for the sake of future reference. 



Lemma 4. Let {Ki, . . . , Kj.} C C[d] he the range of (p^, labeled so that ( |7^ holds, 
and let {Si, . . . , Sr) G V[d] be the corresponding ordered partition, defined by 

{Si,... ,Sr) = {Ki-K2,... ,Kr.l-Kr,Kr). 

Then for each n, 1 < n < r, 

(pZ\Kn) = Sn. 

Proof. We claim that for each i E [d], ipuj{i) is the smallest set Kn containing i, 
a fact which implies the lemma. To prove the claim, suppose i G ^ and 
iPw{i) = Kn, and choose j G [rf] such that (Pw{j) = K^; necessarily j G K^. Then 



so by equation (|T6D, 



forcing wj^" ^ 0. But this implies \Kn\ < \K„\, since by assumption P{Km,j',w) 
holds, forcing Kn = Km- The claim follows. I 
We define vr : ^ V as follows. 

Definition 5. Let w G B^, let {Ki, . . . ,Kr} denote the range of ip^, indexed such 
that ( f/^ holds, and let 

{Si, ...,Sr) = {Ki -K2,..., Kr-l - Kr, K,.) G 

be the corresponding ordered partition. Set 

7i{w) = {w^%es, © ■ ■ ■ © {w^%^sr e lyf 



Note that the inequalities (\i7\} in the definition of B guarantee that 7r{w) is in 
yy(^i,...,Sr) j^g^ ^/{Sl,... ,Sr) _ Furthermorc, because (5*1, . . . , Sr) determines 

the corresponding chain Ki D ■ ■ ■ D Kr via 

Kn = Sn U Sn+l U . . . U Sr, 

the restriction of vr to 7t~^{W^^'"''^''^) is a coordinate projection. Next we define a 
map p : V B^, which will turn out to be the inverse of vr. 

Definition 6. Let .•S'r) ^ {Ki, . . . ,Kr} be the corresponding chain 

in C[d], defined by 

Kn = Sn U Sn+l U ■ ■ • U 5"^ 

for each n, 1 < n < r . Set 
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where w ^ W is defined in terms of its components as follows. Given 5 C [d], 
choose Ki ^ S to be the smallest member of the above chain that includes S , and set 



(19) 



lfn^S^ 



To verify that w = p{w^^^''" '^''^) G B^, we check that equations (|T6|) and inequahties 
( |T7| ) of Definition |^ are satisfied. Indeed w satisfies the inequahties (0) by ( |T9| ) and 

( S S ') 

the definition of ' ■ To see that w satisfies equations (0), let i? C S* C [d] 

with i,j G R. Let ^ Rhe the smahest set in {Ki, . . . , i^^} that includes R, and 
let Kn ^ S he the smallest that includes S; thus i^m ^ Kn- Applying Definition 
the left-hand side of (|T^ works out to be 



(20) 



C'-Vw' = I Ct''iCf-~^wf-){C^--'wh if J e Sm and ^ e S„ 
* * 1^ otherwise 

-iKm—RryKn~R Sm„,,Sn 



Cj^'-'^Cf ""^wj^wf" if J e and t E Sn 
otherwise 



the latter equality since Cf" "^Cf"^ = Cf " ^. The right-hand side of ([16]) is 



C''~Vw' = i C'/"«(Cf ™" V")(C^f ""^^f ) if ^ e 5^ and j G 5„ 
* 1 otherwise 



(21) 



Cfr^~RcKm-R^s^^s^ iii^Sm and j G 5„ 
otherwise 



since Cf ""^C/"^ = Cf"'^. Now, if m ^ n, then 5„ C - /sT^, and therefore, 
since i, j E R K^, neither i nor j belongs to Sn and consequently both ( PO] ) and 
( |?I| ) are 0. On the other hand, if m = n, then (^) and ( |?I| ) are precisely the same. 
So equations (|T6|) are satisfied in any case, confirming that p(^w^^^''"'^''^) G B^. 

Proposition 7. T/ie maj»s p and ir are inverses, that is, p on = Iqa and nop = 1— a. 

Proof [n o p = l:pA.) Let y;('S'iv -Sr) ^ g^,^^ ^ ^ p(ty('^i'--- ••^'■)); we need to 
check that n{w) = w^^^''"'^^\ It follows directly from (|l9|) that for n G S*!, 

V9^j(ra) = /Cj = Si U Si+i U . . . U 5^. 

Then, by Definition ^, 

7r(t5) = {W^')nes^ © ■ ■ ■ © (^f)nG5., 

where, again by (p^ , tl;^' = C^'~^'w^' = so that 7r({i;) = y;('S'iv ,s'r)_ 

{p o n = IgA.) Let w E B^ and let ^/(•S'lv .^r) _ 7]-(^|j}) its image by vr, so that 
= , where Kj = Sj U Sj+i Li. . .USr- In order to verify that p{w^^'^''" ''^'"•') = w, 
fix S* C [d], and, following Definition choose 



i^i = S'i u s'i+i u...uSr^ s 
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to be the smallest such set that includes S, so that S* fl 7^ 0. Then 



(22) 



line Si 
if n ^ 

ifn^^i • 

For comparison, fix m E S H Si, and compute 5;^ using equation ([T6|): 

(23) ^ = (7™?:) -^Sf . 

We will verify that the representations (^) and (^) are proportional, which implies 
that p(^w^^^''"'^''^)^ = , as desired. By Lemma H, if n ^ 5*1, that is, if n G 5*1+1 U 
. . .U Sr, then 

ipain) E {Ki+i, . . . ,Kr} 

and consequently w^^ = 0. On the other hand, if n G Si, then ipw{n) = Ki and w^^ 7^ 
0; in particular, 7^ 0. The coefficient w^/ {C^^~^w^') in (|^) is independent of 
n E S, and non-zero since by definition has at least one non-zero entry, therefore 
(|2^) and ( p3D are proportional. I 



2.4. Stable sequences. We focus our analysis of / : mainly on its values 

at sequences in W^^ , the principal component of T> , that tend to the boundary of 

in Vrl'^l. The blow-up {B^, it) of has a central place in this analysis; indeed, 
as mentioned earlier, we are working toward a proof that 

is a homeomorphism. We keep the notation p = vr"^, consistent with Definition |^ 
and Proposition ^ Let denote the piece of B^ corresponding to W^, that is, 

Bt = p{W^f). 

Note that B^ is characterized as the collection of points w E B^ for which each entry 
w^, n E S C [d\, is non-zero. For w E B^ the range of ipw is the singleton set {[d\}. 
The restriction of tt to B^ coincides with projection of w onto the coordinates w'^', 
and is a homeomorphism with respect to the natural topology of w!f^ . 

Definition 8. Call a sequence w{n) E w]f^ , n E 'L^ , stable if the corresponding 
sequence p{w{n)) E B^ is convergent in W (equivalently, in B^). 

By Definition convergence of p{w{n)) in W amounts to the simultaneous con- 
vergence of the sequences 



d'^~^w(nl) in W^, 

i&S 
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for each non-empty S C [d]. Since the constants Cl'^ ^ {i & S) are strictly positive, 
this is equivalent to the simultaneous convergence of 

{w{n)i)i^s in W^. 

Thus we have the following characterization. A sequence w{n) G is stable if and 
only if, for every i,j G [d], there exists Lij G R"*" U {0, 00} such that 

w(n)i 



w[n)j 



Lji as n — s> 00. 



The Lij determine an ordered partition (5*1, . . . , Sr) G V[d], where cells of the partition 
are defined by the equivalence relation 

i j 4^ < Lij < 00, 

and the ordering of the cells Sm < Sn is determined by 

i E Sm, j & Sn ^ Lij = 00 (equivalently Lji = 0). 

We call {Si, . . . , 5*^) the limiting partition of the stable sequence w{n) G W^l^; it is 
also characterized by the fact that 

TT f lim p{w{n))) G 

\n— »oo / 

Let us use the notation 

^ w^' = w-^i © ■ ■ ■ © 

{Sl,...,Sr) 

for w^'^ © ■ ■ ■ © w^"" G Given a stable sequence w{n) in with limiting 

distribution {Si, . . . , Sr), we calculate 

vr ( lim p{w{n)) ) 

as follows. As usual, for each i, 1 < i < r, set 

Ki = SiU Si+i Li . . .Sr. 

By Definition ^ the restriction of tt to p '^''^^ coincides with a coordinate 

projection on W, call it 7r('^i' - ''^'-). Then 

vr (lim p{w{n))) = -jii^^'-'^^) ( lim p{w{n))) 

\n— ♦00 / \n—>oo / 

= lim ■n^^'-'--'^''\p{w{n))) 

n— >oo 

{Sl,...,Sr) 



= hm m fcfl-^'t^(n),) 

{Sl,...,Sr) 

(24) = lim (Cfl-^^t.(n),) . 



(S'l,...,5r) 
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Returning to distributional notation, write 

dai{n) = w{n)j 5{\ — Xj), 

so that supp dai{n) is indexed by Si, and let 

denote the characteristic polynomial of f{dai{n)), which does not depend on n. The 
assumption that w{n) is stable guarantees that each 

dai = lim dai{n) 

n— >oo 

exists and has support indexed by Si. Rewriting (|2^), we have 

(25) n ( lim p{dai{n) + . . . + dar{n))] = dai ® p^^^ da2 © ■ ■ ■ ©pLi4- +da _i d^r. 
Now, the sequence 

(26) da{n) = dai + da2 + . . . + n~^^^ dcir 

is stable, has limiting distribution (5*1, . . . , 5*^), and p{da{n)) converges to the image 
of ( ^51) by p. Replacing in (|26|) with the continuous variable t, < t < 1, yields 
the following. 

Proposition 9. For every stable sequence w{n) in VF^' there exists a corresponding 
moment curve in of the form 

da(t) = dai + t da2 + . . . + t^~^ dar, <t < 1, 

where {suppdai, . . . , suppdar) G Va, such that 

TT ( lim p{w{n)) ) = TT I lim p{da{t)) ] 

Vn— >oo / \ t^o+ J 

= dai © pLi da2®---® pLi+...+da._i dar- 

Thus in order to study the way the pieces of are embedded in by p, it is 
enough to study the behaviour of moment curves. In particular. 

Proposition 10. = clB^. 

Proof B^ is closed, so cli3^ C B^. To see that the reverse inclusion holds, let 
w & B^ and write 

dai © ■ ■ ■ © dar = tt{w). 
Then w is a limit point of B^, as we can see by choosing the appropriate moment 



curve: 



^ = }}^+ P (^"1 + ^Pdl da2 + ... + f ^Pdl,+...+da^_dar^ 
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Essentially the same argument as in the above proof, carried out within the lower- 
dimensional pieces of V , reveals that has a natural "facial structure" 
identical with that of a special convex polyhedron known as a permutahedron. The 
details of this combinatorial analysis appear in Section |^. 

2.5. Two technical lemmas. In this section we derive two technical lemmas con- 
cerning stable sequences in w]f^ by directly analyzing the maps 



of Definition |l[ 

To begin, we recall some necessary facts and set up notation. Given a matrix 
J, let x{J) denote its characteristic polynomial, and let Jj denote the leading i x i 
submatrix. Note that for J = f{w), w G the sequence {po, . . . ,pd) = "^{w) 

satisfies Pi = x{Ji)^ ^ ^ ^ d. Given w G and a non-empty set Q C [d], let w\q 
denote the restriction of w : [ci] — > R to Q. Recall that, with the notation we have 
established earlier. 



The following notion will be used only fleetingly in the present section; in Sec- 
tion 1^, however, it will be used extensively. A refinement of an ordered partition 
(5*1, . . . , Sr) G is an ordered partition of the form 



to indicate that (-Ri, . . . , Rq) is a refinement of (5*1, . . . , Sr); conversely, we say that 
(5*1, . . . , Sr) is a coarsening of . . . , Rg). 

We note a pair of bounds on the values of on w]^^ , in preparation for the coming 
lemmas. Let P, 5 C R denote the sets 



^ : V 



MOP 





{Ri, . . . , Rg) :< {Si, ... ,Sr) 



P 



{pK^k) I P^ = ^i{w), G 1^1", Afc G A, < z < , 



S 




Proposition 11. Let U = supP and L = inf S. Then U < oo and L > 0. 



Proof. Let S" C R denote the set 



S' 




18 



PETER C. GIBSON 



and observe that 5" ^ S. Now, in S', each ^^ggpf (Afc) > 0, since pi has degree i and 
\Q\ > i. Moreover, the map x{J) is continuous in the entries of J, and the domain 
J' is compact, unhke Therefore 

L' = MS' > 0, 

which imphes positivity of L > L', since S C S'. That U < oo also follows from 
compactness of J7 , by a similar argument. I 
Note that if (5*1, . . . , Sj.) G V[d] is the limiting partition of a stable sequence w{n) G 

w!l^, then for any coarsening {Q, R) G V[d] of (S'l, . . . , Sr), 

lim = WieQ, jeR. 

n^oo w[n)i 

Of course, r > 2 is necessary for such a coarsening to exist. 

Lemma 12. Let w{n) he a stable sequence in with limiting partition {Si, . . . , Sr), 
where r > 2, and let {Q, R) G V[d] be a coarsening of {Si, . . . , Sj.)- If f{'w{n)\Q) con- 
verges, then, for each N , 1 < N < 2\Q\ — 1, the sequence fN{'w{n)) converges, and 



(27) 



lim fN{w{n)) = lim fN{w{n)\Q). 



Proof. We proceed by induction on A^, in the range 1 < N < 2\Q\ — 1. We first 
consider the case N = 1. By Definition |l], 

A + B A 
fi{w{n)) = ^ ^ ^ , fi{w{n)\Q) 



where 



C + D 

A = Y,\w{n)„ B = Y,\ 



Win) 



ieQ 



ieR 



C = ^w{n)i, D = ^ w(n)i. 



In order to show that 

it suffices to show that 
(28) 



lim fi{w{n)) = lim fi{w{n)\Q) 



Ah 


hB 


A 




B 








VD 


~ C 




C + D 


C 


\c + d) 



0, 



as n ^ oo. Now, A/C = fi{w{n)\Q) converges in R by assumption; and we know 
that for every i E Q, j E R, \\mn^oow{n)j/w{n)i = 0, which implies 



lim 



B 



lim 



D 



n—>oo C -\- D n— >oo C + Z) 



0, 



and hence 



SPECTRAL DISTRIBUTIONS AND ISOSPECTRAL SETS OF TRIDIAGONAL MATRICES 19 



Next, suppose that (|27| ) holds for every in the range 1 < < 2z — 3, where 
2 < i < IQI; we will prove that (0) then holds for iV = 2z - 2 and N = 2i - 1 also. 
In the case = 2i — 2, it is simpler to prove the equivalent statement to (0), that 



(29) 
Write 



lim /2i_2(w^H) = lim /2i-2(«^(^)l( 



(p^, ...,p-,) = vl^(^(n)), (go", . . . , gPQi) = vl/(«;(n)|Q), 



where here the superscript n in and is not an exponent, but just an index. By 
Definition |1|, 



A + B 
C + D' 



/2i-2(«^HlQ) 



A^ 



where 



keQ 



k£R 



C = J2(pt2M)Mn)k, D = Y^{pU{\))Mn)k. 

k&Q keR 

^' = Y.^qtM)Mn)k. C = Y.^qt2{^k))Mn)k. 

k£Q keQ 

In order to prove (^), it suffices to show that the right-hand side of the inequality 
(30) 



A + B 


A' 


< 


A + B 


A 




A 


A' 


C + D 




C + D 


" C 


+ 


C 





tends to as n ^ cxD. 

Since p'^_i = x{Ji-i)^ where J = f{w{n)), and g,[!_i = xiJ'i-i): where J' = 
f{w{n)\Q), the induction hypothesis and continuity of x imply that 1 
as n ^ oo. If z = 2, then trivially p"_2 = g"_2; if « > 2, the same reasoning as above 
shows that \C/C'\ ^ 1 as n ^ oo. Since A'/C converges by assumption, it follows 
that 



lim 

n— ►oo 



A A' 

c'c 



and hence that A/C converges. We now argue that the remaining term on the right- 
hand side of (BO), 



(31) 



Ah 


VB 


A 




B A f D \ 




-D 


~ C 




C + D C\C + d) 
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tends to as n ^ oo. For each n, let w{n)^[^ = w{n)^n), where the index k{n) E Q 
renders \w{n)k{n)\ minimal over Q for the given n. Then, by Proposition |Tl|, 

E/w(?2)niin < U^w{n)k/w{n)m:m -> asn^cx), 
L'/u7(n)inin < U^w{n)k/w{n)m:m ^0 as n oo, 

k&R 

C/w{ Ti)ijim ^ L > independently of n. 

It follows that (pTD tends to as n — > oo, proving (|29|) . We omit the proof that 

lim /2i-i(u'(n)) = lim /2i-i(tf («)|q), 

which goes along exactly the same lines as the proof of (^9]), with obvious minor 
modifications. I 

Lemma 13. Letw{n) be a stable sequence in W^, with limiting partition (S*!, . . . , Sr), 
where r > 2. Then 

lim f2\s^\iw{n)) = 0. 

n— »oo 

Proof. To simplify the notation, write = l^il, and set Q = Si and R = S2l-i. . .USr- 
As before, let w{n)^in = w{n)i:(n), where k{n) G Q renders |u^(n)inm| minimal over 
Q for the given n. The fact that 5*1 is a cell in the limiting partition of w{n) implies 
that 

(32) VA; G Q, the sequence w{n)k/w{n)^i^ is bounded 

(indeed it is convergent). Set p'% = \E'Ar(ty(n)) and q"^ = \&Ar(w(n)|Q), where the 
superscript n in and g]^ is an index and not an exponent, so that p^ = x{Jn), 
where J = f{w{n)), and q^ = x{J')i where J' = f{w{n)\Q). In fact q"^ has constant 
value 

independent of n. The entries of Jn, J' are, respectively, 

(/l(^^;(n)), . . . , f2N-i{w{n))), {fi{w{n)\Q), . . . , f2N-i{w{n)\Q)). 
Therefore, by Lemma |12|, lim^^oo Jn = lim„^oo J', and by continuity of x, 

limp-(A)= lim q-{\) = T[{X-\k). 

n— >oo n— too 

Using (p^), this implies 

(33) lim V(p^(Afc))2«;(n)fc/«;(n)niin = 0. 

n— >oo ' ^ 

kdQ 
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For every k E R, w{n)k/w{n)^^i^ — > 0, so, by Proposition [11], 
(34) ^(|?^(Afc))^U7(n)fc/ty(r2)nim < U '^w{n)k/wmm 0, 



k€R keR. 



as 77, ^ oo. On the other hand, since — 1 < IQI, we have, by Proposition |TT 

(35) 5Z(p:^_i(Afc))'u;(n)fc/u;(n)^in > L > 0, 

fceQ 

independently of n, where p^y-i = N-iiw^n)). By Definition ^ 

2l^VN-i{^k))M'^)k + 2^(p^-i(Afc))'w(n)fc 

It follows from (H), (H and (||) that lim^^oo fmi^in)) = 0. I 

2.6. The flip transpose. In this section we explain the origin of the coefficients 
Cf~^ appearing in the equations for in Definition The underlying technical 
fact is a formula involving the fiip transpose. 

The fiip transpose of a (i x (i matrix M is defined to be the dxd matrix having 
entries 

M^{i,j) = M{d+l-j,d+l-i). 

In terms of our convention ([1|) for labeling the entries (ai, . . . ,a2d-i) of a Jacobi 
matrix J, the entries of are simply the reverse sequence (02^-1,... ,ai)- The 
fiip transpose of a matrix induces, in the obvious way, a corresponding operation on 
probability distributions in T>, namely, define the fiip transpose of da G V, written 
da^, to be the spectral distribution of 

ifida)^. 

Note that fiip transposition leaves the spectrum of a matrix unchanged, so that 
supp da^ = supp da. As an operation on distributions, therefore, fiip transposition 
carries into itself; we write for the fiip transpose of w G W^. It follows 
immediately from the definition that for each i, 1 < i < 2d — 1, 

(36) f,{w) = hd-riw""). 

g 

Define constants in terms of the d-element set A = {Ai, . . . , A^} as follows. Set 

p(A) = n(A-A„), P(A) = ^P(A), 

n=l 

and for each non-empty S* C [d], i E [d], set 

Then, hi i e RC S C[d],we have the relation cf /cf = Cf-^. 
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The following result is a translation of ||^, Proposition 3.3] into our current notation. 
Proposition 14. For w G w!f\ is given by the formula 

Wn = — , Vn e [d]. 

This has an immediate implication for stable sequences. 

Corollary 15. Let w{n) be a stable sequence in with limiting partition {Si, . . . , Sr) ■ 
Then the sequence w{n)^ G IV^' is stable and has the reverse- ordered limiting parti- 
tion {Sr, . . . , 5*1). 

Let QU Rhe a partition of [d] , and let w G w!^'^ . Repeated application of Propo- 
sition |14| yields 



\r) 



Ad] 



Cj W j 

'-^i / ieR 

(37) = {C^).,n. 

The following technical result accounts for the coefficients that appear in the equa- 
tions defining B^. 

Lemma 16. Letw{n) be a stable sequence in w!^^ with limiting partition {Si, . . . , Sr), 
where r > 2, and let {Q,R) G V[d\ be a coarsening of {Si, . . . ,Sr)- If f{w{n)^\R) 
converges, then for each i, 2|Q| + 1 < i < 2d — 1, fi{w{n)) converges and 



lim fi{w{n)) = lim fi-2\Q\ ({i^{n)^\R)^) 

n^oo n— »oo 

lim fi-2\Q\ ({Cfw{n)j)j^n 



n— >oo 



Proof. By Corollary 15 w{n) is stable with limiting partition {Sr, ■ ■ ■ ,Si), of 
which {R,Q) is a coarsening. If f{w{n)^\fi) converges then, by Lemma 0, for each 
^, 1 < ? < 2|_R| — 1, fi{w{n)^) converges, and 

lim fi{w{nY) = lim fi{w{n)^\R). 

n— too n— ►oo 

Applying (|36|) to the left- and right-hand sides yields 

lim fj{w{n)) = lim fj_2\Q\{{w{n)^\R)^), 

n— >oo n—>oo 

where j = 2d — i and 2|Q| + l<j<2(i— 1. The calculation (|T7D then gives the last 
expression in the conclusion of the lemma. I 
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Lemma 17. Letw{n) he a stable sequence in W^"^ with limiting partition (5*1, . . . , Sr), 
where r > 2. Then 

lim /2d-2|5.|(wH) = 0. 

n— »oo 

Proof. Applying Lemma Oto the sequence w{n)^ yields 



lim/2|5.|(«^H^)=0. 

n— >oo 

By ( ^61) this is equivalent to the conclusion of the present lemma. 
2.7. The main result. In this section we prove that 



is a homeomorphism. 

The proof involves stable sequences in , where Q ranges over non-empty subsets 
of [d]. We use the characterization of stable sequences in terms of limits, to avoid 
dealing with semi-algebraic sets of the form where F C A. Thus a sequence 
w{n) G is stable if, for every i,j G Q, there exists Lij G R"*" U {0, oo} such that 

w(n)i 

L,,- as n — s> oo. 



w{n)j 



Note that if w{n) is a stable sequence in Wj^ then for Q ^ [d\, the restriction w{n)\Q 
is automatically a stable sequence in . 

Proposition 18. Let w{n) he a stable sequence in W^"^ with liraiting partition (5*1, . . . , 5*^,) 
Then f{w{n)) converges, and 

(38) lim fiwin)) =J I ffi lim (C-'''-'''^'''win) 

\{Sl,...,Sr) 

Proof. We argue by induction on = \S\, where S C [d]. If 1 5*1 = 1 then the 
proposition, with Wf in place of w]^, holds trivially. 

Suppose that 2 < N < d, and that for every non-empty Q C [d] of size \Q\ < N, 
the proposition holds with in place of w]^. We will show that the proposition 
then holds for Wf whenever S C [d] and 15*1 = A^. Thus let S C [d] have size IS"! = A^, 
and let w{n) be a stable sequence in Wf with limiting partition (5*1, . . . , Sr) G Vs. 
If r = 1 then w{n) necessarily converges in since / is continuous on (with 
respect to its natural topology), this implies 

lim f{w{n)) = f ( lim w{n) ) , 

n— »oo \n— >oo / 

which conforms to (P^. Next consider the non-trivial case r > 2. Set Q = Si U 
. . . U Sr-i and R = Sr. Then w{n)\Q is a stable sequence with limiting partition 
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(5*1, . . . , Sr-i), and by the induction hypothesis 

hm /(^HIq) = 7 I hm (Cf ^-^^-^^(n) 

n— >oo \ ^-J — ?^— +00 \ 

In particular, f{w{n)\Q) converges, so by Lemma |12] (with S in place of [d]), for every 
M in the range 1 < M < 2\Q\ — 1, fM{w{n)) converges, and 

(39) lim fM{w{n)) = lim fM{w{n)\Q). 

n— »oo n— >oo 

Corollary |l^ implies that w{n)^\R converges in W^, since i? = S'r is a cell of the 
limiting partition of w{n). Therefore f{w{n)^\fi) converges, and by Lemma |16|, for 
M in the range 2\Q\ + 1 < M < 2\S\ — 1, fM{w{n)) converges, and 

lim fMiwin)) = lim /m-2|s| ({C?w{n)i)i^Sr) 

n— »oo n— >oo \ / 

(40) = fM-2\s\(lnn{CS{n)i)^ 



Vn— »oo 



the latter equality since / is continuous on W^. The behaviour of f]\i{w{n)) for the 
remaining index M = 215*1 — 2|S'r| is determined by Lemma 

(41) lim /2|sh2|5.|(w^(^^)) = 0. 

n— »oo 

The statements (^), (0) and (^T]) collectively imply that lim/(w(n)) has the form 

,,SiU...us'j_i^ 



lim f{w{n)) = f\ lim (cf 
/ f lim (C^w{r 



7 ( © 



ieSr 

SiU...uSj-i 



Win) 



Thus the proposition holds for W^, completing the induction. 
Corollary 19. Let w{n) he a sequence in that converges in B^. Then 
(42) lim / o 7i(w(n)) = / o vr ( lim 



Vn— »oo 



Proof. Write w{n) = p{w{n)), where w{n) = 7i{w{n)) is by definition a stable 
sequence in wj^^, and let (5*1, . . . , Sr) denote the limiting partition of w{n). Now, 



f o n{w{n)) = f{w{n)), and, by Proposition |T^, the left-hand side of (|4^) is 

,,SiU...USj_i^ 



(43) lim f{w{n)) = f\ lim (q 

\{Sl,...,Sr) 



Win) 
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On the other hand, by the calculation (p4D, 

71 ( lim {^{n)) = ff^ lim (c^^^"'^^^'^w{n)i) , 

{Sl,...,Sr) ' 

SO the right-hand sides of (^2|) and (|^) are equal. I 
It is a small step from the above corollary to our main result. 

Theorem 20. The map / o tt : is a homeomorphism. 

Proof. For brevity write g = f o ir. Note first that g is a. bijection since each of 
its factors / and vr is. To see that g is continuous, let w{n) G be a convergent 
sequence. It follows from Corollary [l^ and the fact that B^ = clB^ (Proposition |10D 
that there exists a convergent sequence w'{n) G B^ such that 

(44) lim w'{n) = lim w{n) and lim \\g{w'{n)) — g[w{n))\ \ = 0. 

n—*oo n— ►oo n— ►oo 

By Corollary |19|, 

lim g{w'{n)) = g ( lim w'{n) ) , 

n— >oo \n— >oo / 

which, by (|i^), implies 

lim g{w{n)) = g ( lim 



\n— ►oo 

Thus g is a continuous bijection; since B^ is compact it follows that g is a homeo- 
morphism. I 

2.8. The solution to Problem 1. With Theorem ^ we have a complete solution 
to Problem 1. Since the composition 

is a homeomorphism, the topology that J'^ induces on via / is precisely the 

topology that B^ induces on V via vr. To put it succinctly, the blow-up {B^, vr) 

itself constitutes a solution to Problem 1. A sequence w{n) G V converges in the 
induced topology if and only if 7i~^{w{n)) converges in B^, and in this case the 
limiting value may be computed as 

(45) TT ( lim 7r^^(w(n)) ) , 

\n— >oo / 

which does not require evaluation of /. The point is that vr and p = vr"^, unlike /, 
are easy to calculate by hand. 

More explicitly, a sequence w (n) in the principal component w!^^ of V is conver- 
gent if and only if it is stable. (The characterization of convergent sequences in the 

other, lower-dimensional pieces of V , 



w^^{n) G W, 

{Sl,...,Sr) 



{Sl,...,Sr) 
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is virtually the same: each w^^{n) must be stable in W^\) Recall that the essence 

of the computation (^5[) is captured by the limiting values of moment curves in 
as in Proposition In distributional notation, given dai^ . . . , dur G T) such that 

(supp dai, . . . , supp dar) G Va, 

the moment curve da(t) = dai + t da2 + ■ ■ ■ + dar tends to the point 



TT ( hm TT-\da{t)) ] = f ^ (lim f{da{t)) 



= dai © pa^^ da2®---® Pdai+...+dar-i dar. 

The example discussed in Section is a simple instance of this general limiting 
behaviour. 



3. ISOSPECTRAL SETS OF TRIDIAGONAL MATRICES 

We turn now to a particular application of the material in Section ^ deducing the 
topology of the set of symmetric, tridiagonal matrices having prescribed spectrum 
A. This results in some rather interesting discrete geometry. As mentioned in the 
Introduction, the isospectral sets we consider have been studied before — our purpose 
is to provide a new approach. 

We state precisely the problem of interest in Section Then in Section p.2| we 
analyze the combinatorial structure of the semialgebraic set B^, constructed earlier 
in Section ^ We show that has a natural "facial" structure, corresponding to the 
individual pieces of V , and that the facial lattice of B^ is represented by the partial 
order {Va, The combinatorial structure of B^ allows one to very easily construct 
a polyhedral complex, which we label P^, over the manifold of symmetric, tridiagonal 
matrices that have spectrum A. We describe this construction in Sections |3.3| and 
|3.4] . Finally in Sectio n^75| we carry out some explicit computations. We explore the 
particular example of P3: topologically it is a two- holed torus, and its Petrie polygon 
uses every edge. Concerning the sequence of complexes P^, we prove the remarkable 
fact that hyperbolic tangent generates the Euler characteristics. 

We follow the nomenclature for convex polytopes given in Br0ndsted's introductory 
text M. In particular we use the characterization of a proper face of a d-polytope P 
in R as the intersection of P with a supporting hyperplane in R'^. 
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3.1. Symmetric, tridiagonal matrices; Problem 2. Let T denote the set of 
finite-dimensional, real, symmetric, tridiagonal matrices. We consider the following 
problem. 



Problem 2 Let A = {Ai, . . . , A^} C R be a rf-element set. Describe the topology 
of the isospectral set 

r^ = {T eT I specT = A}. 



The set C consists of those matrices whose next-to-diagonal entries are 



non-negative, and so the parametrization 



A 



established in Section ^ applies directly to Problem 2. Indeed let £ denote the set of 
(i- vectors of the form 

e = {ei, ... ,ed), where each ej G {+1, — !}• 

For each e E S let denote the operation ondxd matrices of conjugation by diag(£), 
so that 

Se{M) = diag(£)Mdiag(£). 

Then each map of the form 

is a homeomorphism from onto a "patch" of T^, and collectively the maps Sgo/ovr 
cover T^. (Note that the are a way of encoding the modification of Definition |l| 
which would correspond to a particular choice of positive or negative square root for 
each off-diagonal component of / — as Definition |l| stands, only the positive square 
root is used.) The topology of is thus captured by the complex {B^ x ^)/ ~5 
consisting of copies B^ x {e} of B^, glued together according to the relation 

{w, e) ~ (w', e') <^=^ o / o Tr{w) = v ° / ° '!t{w'). 

This gluing map respects the inherent "facial" structure of B^, which we describe in 
Section below. In principle, the complex {B^ ^ £)/ ~ is all we need in order to 



calculate the standard topological invariants of T . However, to make the combina- 
torial structure of {B^ x ~ more transparent, we go a step further in Section 3^ 



and replace B'^ with a particular convex polytope, call a permutahedron, which is 
combinatorially equivalent to B^ but nevertheless a simpler object. The final result 
is a polyhedral complex over which not only encodes the topology of T^, but 
which is itself rather interesting. 

3.2. The facial structure of B^. For each ordered partition (S*!,... , 5*,) G 
set 
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We refer to B^^^ ^ ^ as a face of B^, in analogy with the terminology for convex 
polytopes. We give partial justification for this later in the present section, by showing 
that, in the appropriate affine setting, each B^^^ ^ is the intersection of B^ with 
the bounding hyperplane of a supporting half-space of B^, exactly as for a convex 
polytope. But B^ is of course not a polytope; it's faces are in general curved. In 
the next section we show that the "curved polytope" B^ can be flattened, that is, 
mapped to a bona fide convex polytope that has the same facial structure. Our main 
objective in this section is to ascertain the inclusion relation between faces of B^, 
and, more precisely, to show that (V[^,^) represents the facial lattice of B^. This 



entails what is in essence a technically messier version of the proof of Proposition ^ 
(which states that = B^). 

To begin, we note the following extension of Proposition ^ 

Proposition 21. Let S C [d], {Si, . . . , S",) G Vs, and w^'^''- G ^yj^i'-'^'-). Then 
there exists a moment curve w{t) e Wf such that 



r' (^\unf{w{t))^ 



Proof. By Theorem the topology that J'^ induces on via / agrees with the 
topology induced on V by B^ via vr. In the case S = [d], Propositions |^ and p!0| 
therefore imply the present proposition directly. But A, and in particular d = |A|, is 
arbitrary, so once the proposition holds for S = [d], it holds for all S G [d], since we 
may replace A with A' = {A„ G A | G S}. I 
For each (5*1, . . . , Sr) G V[d], let -F(Si,...,s^) C B^ denote the set 



,(i?i,...,R,)r<(5i,...,5.) 



(-Rl,... ,Rq) 



Lemma 22. For each {Si, . . . , 5,.) G V[d], ,5,.) ^ ^{Si,...,Sr)- 

Proof. Let w G -F(5i,... ,5^), and write w^^^'---'^''^ = n{w), so that {Ri, . . . ,Rg) ^ 
(5*1, .. . , Sr) and hence has the form 

{Rl, . . . , Rq) = {Rl, ... ,Rl^,... ,R\,... , R^J, 

where each {R{, . . . , R?^^) G Vs^- Using Proposition for each Sj, let w^^{t) G W^^ 
be a moment curve such that 

(\imf{w'^^{t))\ =w^'^'-''^^^\ 
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where w^^^''"'^"^^ 



denotes the appropriate component of w^^^'^^^ '^''^ Then 0(5^ Srf^^^ (^) 



is a moment curve in W_ 



{Sl,...,Sr) 



and 



whereby 



hm f 

t-,o+ 



w'^it) I = ^hm fiw'^it)) 



ASl,... ,Sr) 



{Sl,...,Sr) 



\unf{w'^it)) 



iSl,...,Sr) 



(Sl,...,5r-) 



iRl,...,Rq) 



,{Sl,...,Sr) 



By Theorem 2C, this imphes 



nmvr-i ^^.(t) 

(51,... ,5^ 



SO that w e 'B(\.^„.,,s.) = cl [tt-^ (w^f''-'^''^ 

In Section ^ it was convenient to have in projective space so that we could 
use homogeneous coordinates, and not have to worry about rescahng. For present 
considerations, it is better to embed in affine space. To that end, for each non- 
empty S C [(i], let 

= {w^ -.S ^ R} 

be the |S'| -dimensional Euclidean space of real- valued sequences indexed over S, and 
set 

E= H R'. 

<D^SC[d] 

Let A denote the affine subspace of E defined as 



A 



w e E 



wz 



Now, the ambient space ofB^, W = n0^5c[d] is a quotient space of £' — {0}, and 
the inequalities ( [T7| ) defining B^ guarantee that each w G B^ has a unique normalized 
representor w ^ A. Hence we may regard B^ G A as sitting in affine space, without 
altering its topological structure; we do so for the remainder of the present section. 
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We now define a family of half-spaces H'^^^ of A, along with their bounding 
hyperplanes H(^Si,...,Sr)- 1^ follows from the the explicit construction of p = vr"^ given 
in Definition | that for w G tt~^ (w!^''-'^''^^, the set 

1= {{n,S) \ neSC [d], w^ = 0} 
does not depend on the particular point w. Rather, 

^ = ^{Sl,...,Sr) 

is determined by the partition (5*1, . . . , Sr) G Moreover, 

(46) IiSu...,Sr)^I(Ru-,R,) ^ {Ru...,Rg) ^ {Su...,Sr). 

Now, given {Si,... , Sr) E V[d], let 

H(Sl,...,Sr) C A 

denote the hyperplane in A determined by the equation 

(n,5)Gl(Sj,...,s^) 

and let H^g_^ C A denote the corresponding half-space, determined by the in- 
equality 



Lemma 23. For every {Si, . . . , Sr) G V[d], the half-space Hf^g^ ^^-^ includes B^, and 

-^(5i,...,5.) = H(^Sl,...,Sr) (^B^- 

Proof. The inequalities (|1^ defining and the defining equations of A imply 
that if w E B^, then its normalized representor w E A satisfies > for every 
nE S C[d]. It follows that B^ C H^g^ 

Also, w E H(^Si,...,Sr) if only if for each {n, S) E X(^Si,...,Sr): — 0- (H) ^^i^ 
condition is equivalent to the condition that 

Ti{w) E 

for some {Ri, . . . , Rq) < {Si, . . . , Sr), or in other words, w E F(^Si,...,Sr)- ' 
Combining Lemmas ^ and ^ yields the following characterization of the faces 

Proposition 24. For every {Si,... , Sr) E V[d], B^^^^^^g^^ = F(^Si,...,Sr)- 

Proof By definition, F(^Si,... ,Sr) ^ ,Sr)^ ^ Using Lemma ^ F'{Si,...,Sr) = 

B^ n H(^Si,...,Sr) is closed, so 

Lemma 52] gives the reverse inclusion. I 
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Corollary 25. Sf^,^,... C iSf^^^ , z/ and only if {R,, . . . , Rg) ^ {S,, . . . , S,.). 



Thus we have shown that {V[d], ^) represents the facial lattice of B . Note that the 
dimension of the face Bjg^ ^ ^ is dim{W^^'"''^''^) = d — r. The least-dimensional 
faces, corresponding to ordered partitions of the form 

(Mi)},...,Mrf)}), 

where a is a permutation of [d], consist of a single point and have dimension 0. If 
two ordered partitions 

(i?!, . . . , Rq) :< {Si, . . . , 5*,.) 
are adjacent with respect to ^, that is if r = g + 1, then evidently 



dim(Sf<,^,„„5^))= dim(i3f^^, 
a fact which we will use in the next section. 



+ 1, 



3.3. Permutahedra. Permutahedra are convex polytopes that have a very special 
structure. For instance, they are simple zonotopes that tile Euclidean space; see 
Ziegler [0, pp. 17,18,200]. We give a slightly different description from that of Ziegler, 
and consequently our characterization of faces. Proposition |2^, differs somewhat from 
his. But it is easy to reconcile the two descriptions, so we do not bother to prove 
Proposition 

We fix notation as follows. Let S[d] denote the symmetric group, consisting of all 
permutations of [d] = {1, . . . , d}. For each non-empty subset i? C [rf], let Sr denote 
the subgroup of consisting of those permutations which fix each element of [d] — R; 
let Q denote the collection of subgroups of S[cj] given by 

g = {Sr^ x...xSr^ I ({1}, . . . , {d}) ^ ...,Rq)e } . 

We say that a non-empty set R C [d] is an invariant set of a group H C S[d] if every 
a & H fixes R setwise; if in addition no non-empty proper subset of R is invariant, 
we say i? is a minimal invariant set of H. Note that for 

G = SR,x---xSR^eg, 

the sets Ri, . . . ,Rq are the minimal invariant subsets of G, and they have a unique 
ordering for which ({1}, . . . , {d}) ^ . . . ,Rg) G V[d]- The collection Q will be 
useful in discussing the structure of the permutahedron, which we now define. 

Definition 26. Let Pd denote the convex hull of the set 

(47) {(a(l),... ,a(rf))GR'^ | a E S^d]} ■ 

The convex polytope Pd is called a permutahedron. 

Given a set F C R'^, let F denote its convex hull. We will generally identify 
elements of S[d] with points in R"^ as in (^7\); with this convention we have, for 
instance. 
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The set of groups Q provides a simple description of the faces of Pd-, as follows. Let 
QS\^d] denote the set of right cosets of the form Gcr, where G E Q and a G S\^d]- 

Proposition 27. The faces of Pd are precisely the sets F where F E GS\^d]- 
The above proposition implies that 

i^Sid],^) 

represents the facial lattice of Pd- Indeed this corresponds to the facial lattice of 
via the following map. 

Definition 28. We denote by if the map 

^ ■ GS[d\ V[d\, 

constructed as follows. Let H G QS[d]. Choose a E H arbitrarily; then Ha^^ G Q. 
Let Ri, . . . ,Rq C [d] be the minimal invariant sets of Ha"^ , ordered such that 

({!},..., {4) 

Set 

Note that if if G QS[d] corresponds to a vertex of Pd, that is, if if = {a} is a 
singleton, then according to the above definition 

^{H) = {{a-\l)},...,{a-\d)}). 

In other words, on the level of vertices, ip coincides essentially with the map cr 
a~^. More generally, let H G QS[d] and write (Ti, . . . ,Tr) = ^{H). Then it is a 
straightforward to verify that 

(48) H={ae S[d] I ({1}, . . . , {4) ^ (^m), • • • , ^(Tr)) } , 

whence ip is one-to-one. And given any (Ti, . . . , T^) G Vid], the set H defined by (^81) 
belongs to GS[d\, and so p is onto; that is, p is a bijection. Indeed, Pd and are 
combinatorially equivalent, in the following sense. 

Proposition 29. The map ip : QS[d] — ^ Vid] is a lattice isomorphism, so that 

(^5[,],C)-(P[,],^). 

Proof. Since p is a. bijection, we need only show that for every H,K E QS[d], H (1 K 
if and only if ^{H) ^ fiK). To that end, let H,K E QS[d] and write (i?i, . . . , Rg) = 
p{H), (Ti, . . . ,Tr) = p{K). By (|48| ) the condition H O K is evidently equivalent to 
the condition that for every a G S[d], 

({1}, . . . , {4) ^ (a(i?i), . . . , a{R,)) ^ ({!},... , {4) ^ {a{T,), a{Tr)). 

And the latter condition is easily shown to be equivalent to the relation 

(i?i, . . . , Rq) :< (Ti, . . . , Tr). 

I 

For an arbitrary convex polytope P, if two faces F d F' are adjacent in the facial 
lattice of P, that is, if there is no face F" of P such that F G F" G F', then 
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dim(F)' = dim(F) + 1. As pointed out in the previous section, the facial lattice of 
also has this property; therefore the map ip automatically respects dimension: 

dim(F) = dim(i3^(^)), 

for every F G GS[j\. 

Any isomorphism l between the facial lattices of two convex polytopes P and P' has 
a topological realization, in the sense that there exists a homeomorphism h : P ^ P' 
which preserves faces, and whose action on faces agrees with l. By the same token, 
the map ip has a topological realization 

h:Pa^ B^. 

We sketch briefly one possible construction of such an h. Here it is useful to have a 
metric on B^, so that we may speak of geodesies. As in the previous section, let us 
regard B^ as sitting in the affine space A, and consider the metric that B^ inherits 
from the ambient Euclidean space E = Yl^^scid] R-''^'- 

For each face B^g_^ ^ ^ of B^, consider the point ^^('S'l'--- .-Sr) ^ yj/jSi,... ,5r) j^g^yj^g 

constant entries Wn = 1 for every j, n E Sj. We refer to the normalized representor 
b{Si,...,Sr) G A of 7r~^(w('^i'--- ''^'■)) as the "barycentre" of B^g^ g^y The map h may be 
constructed recursively as follows. 

Initial step For each vertex a G Pd, set h{a) = -Bj^do-}) (the latter consists of just a 
single point). 

Continuing step Suppose that h has been defined on the j-skeleton of Pd, that is, 
on the set of faces of P^ having dimension at most j, where < j < — 2. Given a 
(j + l)-dimensional face F, extend /i to F as follows. Let h map the barycentre /o 
of F (the mean of F) to the barycentre of ^!^{f)- And for each point fi in the 
j-skeleton of F, let h map the segment /o/i to the geodesic in B^^^ ^ ^ connecting 
bip{F) to h{fi), at constant speed with respect to arclength. 

The map h : P^ ^ B^ so constructed is a homeomorphism, and its inverse can be 
viewed as a "flattening" of B^. The permutahedron Pd is homeomorphic to via 
the composite map 

(49) p^JUB^'^v''^ j"". 

This serves as the basis for a polyhedral complex, which we describe in the next 
section, on the set of d x d symmetric, tridiagonal matrices having spectrum A. 

3.4. A polyhedral complex over T^. Referring to the composite map (^9]) , for 
p E Pd write 

Jp = f OTTO h{p). 

Let S C R'^ denote the 2'^-element set consisting of maps e of the form 

e:[d]^{+l,-l}. 
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(In other words, £ consists of the vertices of the twice-unit cube centred at the origin.) 
For each (p, e) E Pd^ £, set 

g{p,e) = diag(e) Jpdiag(£). 
Since the map p t-^ Jp is a homeomorphism from onto J'^, the map 

g:PdX£^r^ 

is easily seen to be a continuous surjection. Let ~ denote the relation defined on 
PdX £ defined by 



And let [{p, e)] denote the equivalence class of {p, e) modulo ~. The relation ~ may 
be characterized purely in terms oi Pd x £ as follows. 

Proposition 30. Let {p, e), {p', e') E PdX £ . Let F G QS[d] represent the unique face 
F of Pd whose relative interior contains p. Choose a E F, and write ... , Rq for 
the minimal invariant sets of the group Fa~^ G Q , ordered such that 



Then {p,e) ~ {p',£') if arid only if p = p' and £^|_Rj = i^'lfij for each Rj. 

Proof. The matrix Jp may be obtained from g{p, e) by replacing off-diagonal entries 
with their absolute values. Since p ^ Jp is one-to-one, it follows that g{p,e) deter- 
mines p uniquely. Thus (p, e) ~ (p', e') implies p = p' and we need only characterize 
the relation {p,e) ~ (p, 

The fact that p is in the relative interior of F implies that h{p) is in the relative 
interior of ^^^i^py or in other words, in vr"^ [w"^^^^^. Therefore Jp E f {w'f"^^^, 
which, recalling that '^{F) = ((7~^(i?i), . . . ,a^^{Rq)), implies 

t/p = t/i © ■ ■ ■ © t/g, 

where each Jj is an x Jacobi matrix. With this decomposition, the condition 
(p, e) ~ (p, e') becomes 

diag(£:) Ji © ■ ■ ■ © Jg diag(£) = diag(£') Ji © ■ ■ ■ © Jg diag(£:'), 

which is evidently equivalent to the condition that, for each Ri, = ±e'\R-. I 
We now define the polyhedral complex of interest. 

Definition 31. Let Pd denote the quotient {Pd x £)/ ~, and define 



to be the map induced by g, so that g ([(p, e)]) = g{p,e). 

Faces of Pd are sets of the form F x {e} where e E £ and F is a face of Pd- 
Proposition |30| implies that either the relative interiors of two faces F x {e} and 
F' X {e'} are disjoint, or the two faces are equal. And the latter can only happen 
if F = F'. In the case of a least-dimensional face, a vertex {a} where a E S[d], we 
have {a} x {e} = {a} x {s'} in Pd for every e E £. In other words, Pd has the same 




{{1},... ,{d})^{Ru... ,R,). 



g-Pd 
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vertices as P^, one for each permutation a E S[d]- At the other extreme, there are 
2^^"^ distinct faces of Pd of the form x {e}; a given set Pd x {e} gets identified in 
Pd with the set Pd x {e'} only when e' = ±e. 

The following result encapsulates the solution to Problem 2. To prove it, we need 
only verify that the map g is a. homeomorphism. 

Theorem 32. Pd is a polyhedral complex over . 

Proof. Since g : Pd x S is a continuous surjection, ^ is a continuous bijection 

by construction. That Tj is a. homeomorphism then follows from compactness of Pd- 

_ I 

Thus the topology of is encoded by the polyhedral complex Pd- Using Proposi- 



tion 30, the calculation of topological invariants of Pd is reasonably straightforward; 



we illustrate this in the next section. 



3.5. Illustration of the results, Euler characteristics. In this section we apply 
the earlier results specifically to 3 x 3 matrices, and we calculate the Euler charac- 
teristic of for an arbitrary d-element set A = {Ai, . . . , A^} C R. 

To begin, consider the permutahedron P^, the basic unit of the complex P^, and 
let us count the number of faces having a fixed dimension n, where < n < c? — 1. 
A face F of Pd has dimension n if and only if the ordered partition 

has r = d — n cells. There are precisely 

(50) {d-n)\S{d,d-n) 

such ordered partitions, where S{d, d — n) denotes the Stirling number of the second 
kind. 

Recall that if F,F' G GS[d] are distinct, then for any e,e' E £ so are the faces 
F X {e},F' X {e} of Pd- On the other hand, for fixed P, Proposition ^ tells us 
precisely how many faces of Pd there are of the form P x {e} [e E £)- Namely, let 
N denote the number in question, and write 

{a-\R^),... ,a-\R,))=ip{F), 

where a E F and Pi, ... ,Rq are the minimal invariant subsets of Pcr^^. Then is 
the number of equivalence classes of S modulo the relation 

€\r. = ±€'\r^ for every Rj. 

That is, 

(51) = 2^/2^ = 2", 

where n = d — q is the dimension of P x {s}- Combining (|50D and ( ^I]) yields the 
number of n-dimensional faces of Pd and hence the following explicit form for the 
Euler characteristic. 
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132 123 




Figure 1. The complex P3 with an individual facet depicted sepa- 
rately at left. The letters a-h and their negatives indicate the gluing 
rules for edges of the 16-gon. There are 12 distinct edges (including 
the four interior ones), 6 vertices and 4 faces. 



Theorem 33. The Euler characteristic of Pd is 

d~l d 

x(Pd) = $^(-2)"(d - n)\S{d, d-n) = 5^(-2)'^"«g!5(d, q). 

n=0 q=l 

We show below that the above value of the Euler characteristic has a simple gen- 
erating function. But first we discuss the particular complex P3. It consists of 
2^"^ = 4 copies of the regular hexagon P3, glued together according to the relation 
~ as depicted in Figure |l]. By Theorem ^ the Euler characteristic is 

X (^3) = -2, 

which, using Theorem ^ and the standard classification of 2-manifolds, implies: 

Theorem 34. is a two-holed torus whenever K is a set of 3 distinct real numbers. 

From the point of view of discrete geometry, the polyhedral complexes Pd seem 
to be of intrinsic interest, quite apart from their relationship to tridiagonal matrices. 
For instance, J. Bokowski has made the following observation. See pp. 24, 25, 223] 
for relevant definitions and comparison with the Platonic solids. 

Proposition 35. The Petrie polygon of P^ covers every edge. 
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This unusual property can be verified directly in Figure |I]. Another remarkable 
property of the complexes Pd — and also of the isospectral sets — is that hyperbolic 
tangent generates the Euler characteristics. 



Theorem 36. 

oo 

tanhx = X {Pd) 



d\ 

d=l 



Proof. According to a theorem of Frobenius |]^, p. 244], the so-called Eulerian 
polynomial of degree d has the representation 

d 

pa(u) =uY, r^-S{d, r) (n - 1)^"''. 



r=l 



It is known [|, p. 259] that 

oo 

tanx = ^\Pd{-l)\'^ 



x^ 



d=l 

oo 



x'^ 



d=i 

where i"^ = —1. It follows from the relation 

tanhx = — ztan(zx) 

that 



x"^ 



tanhx = 



d=l 

oo 



d=l 

the latter equality because tanh is odd. Comparing Theorem |3^ to the Frobenius 
representation of pd shows —pd{—l) = x{Pd)- I 
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